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^ . We study the orbital stability and instability of single-spike bound states of semi- 

I classical nonlinear Schrodinger (NLS) equations with critical exponent, linear and non- 

linear optical lattices (OLs). These equations may model two-dimensional Bose-Einstein 
p condensates in linear and nonlinear OLs. When linear OLs are switched off, we derive the 

\ asymptotic expansion formulas and obtain necessary conditions for the orbital stability 

and instability of single-spike bound states, respectively. When linear OLs are turned on, 
we consider three different conditions of linear and nonlinear OLs to develop mathematical 
theorems which are most general on the orbital stability problem. 

1 Introduction 

. Recently, optical lattices have created many interesting phenomena in Bose-Einstein con- 
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densates (BECs) and attracted a great deal of attention. Two types of optical lattices are 
considered: a linear optical lattice (OL) (cf. |28]) and a nonlinear OL (cf. [1] and |35]). A linear 
O I OL is a series of potential wells having a periodic (in space) intensity pattern which may confine 
Q ■ atoms of BECs in the potential minima. A nonlinear OL can be obtained by inducing a periodic 
spatial variation of the atomic scattering length, leading to a periodic space modulation of the 
. nonlinear coefficient in the Gross-Pitaevskii equation (GPE) governing the dynamics of BECs. 
^ ] The GPE is a nonlinear Schrodinger (NLS) equation in the presence of the Kerr nonlinearity 
] describing a BEG in a linear and a nonlinear OL given by 

-^^ = DA,p - Vtrapi^ - g\ij\'ij , (LI) 

for X G M^, N < 3 and t > 0. Here ip = ip{x, t) G C is the wavefunction, D is the diffraction (or 
dispersion) coefficient, and Vtrap is the potential of the linear lattice. Besides, g = fim{x) ~ a 
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characterizes the nonhnear lattice, where a denotes the spatially modulated scattering length, 
/i is a nonzero constant and m(x) = m{xi, ■ ■ ■ ,X]^) > is a function depending on spatial 
variables (transverse coordinates) Xi, ■ ■ ■ ,xj\f (cf. [2], [O])- 

The underlying dynamics of fll.ip is dominated by the interplay between adjacent potential 
wells of linear OLs and nonlinearity of nonlinear OLs. When the nonlinearity is self-focusing 
i.e. D > and /i < 0, a balance between these two effects may resist collapse or decay and 
result in bright solitons. Experimentally, bright solitons can be observed in linear and nonlinear 
OLs, respectively. One may find stable bright solitons in three-dimensional linear OLs (cf. [7]). 
On the other hand, two-dimensional bright solitons can also be investigated in two-dimensional 
nonlinear OLs (cf. [13]). Consequently, under the influence of linear and nonlinear OLs, two- 
dimensional bright solitons must have suitable stability for experimental observations. However, 
most theoretical results (e.g. [10] and [H]) focus on the orbital (dynamical) stability of only one- 
dimensional single-spike bound states which are steady state bright solitons in one-dimensional 
nonlinear OLs without the effect of linear OLs. To see how linear and nonlinear OLs affect the 
stability of two-dimensional single-spike bound states, we develop mathematical theorems for 
the orbital stability and instability of two-dimensional single-spike bound states of (11. ip under 
different conditions of linear and nonlinear OLs. 

To get two-dimensional single-spike bound states of (11. ip . we may assume N = 2, D > 
and the scattering length a, i.e., n is negative and large due to the Feshbach resonance (cf. [l]). 
Setting h"^ = D/{—fi), V{x) = Vtrap{x)/ {—fi) and suitable time scale, the equation (II. ip with 
negative and large /i can be equivalent to a semi-classical nonlinear Schrodinger equation (NLS) 
given by 

-ih^ = h'^A'ilj-V^lj + m\tlj\'^'ilj, xGM^t>0, (1.2) 
dt 

where < /i ^ 1 is a small parameter, V = V{x) is a. smooth nonnegative function and 
m = m(x) is a smooth positive function. For the spatial dimension > 1, we may generalize 
the equation (II. 2p to a NLS having the following form 

-th^ = h'^Aip-Vip + mMP-'^tlj, xGM^,t>0, (1.3) 
at 

with critical exponent 

p = l + l, N>1. (1.4) 

In particular, when N = 2, the equation (II. 3p with (II. 4p is exactly same as (II. 2p . 

Single-spike bound states of (II. 3p are of the form ijj{x,t) = e^^^^'^u{x), where A is a positive 
constant and u = u{x) is a positive solution of the following nonlinear elliptic equation 

h^Au- {V + X)u + muP = 0, ueH\R^), (1.5) 

with zero Dirichlet boundary condition, i.e., u{x) — )■ as \x\ — t- oo. When = and m = 1, 
problem (11.51) admits a unique radially symmetric ground state which is stable for any A > 
if p < 1 + jf, and unstable for any A > if p > 1 ^ (cf. [1], [8] and [13]). For V ^ or 
m ^ 1, there exists Uh a single-spike solution of (II. 5p . provided both V and m are bounded and 
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satisfy another conditions, for example, conditions in the following Theorem ll.mi.4l (cf. [20]). 
For other other nonlinearity in the possibly degenerate setting, see [5], [H], [12], [SI], [32], [57] . 
[5U] . [inj, [H] and reference therein. Hereafter, we set iph{x,t) := e''^^/'^Uh{x) as a single-spike 
bound state of (11.31) . where Uh is the single-spike solution of (II. 5p . 

In this paper, we want to study the orbital stability of the bound state iph for the equa- 
tion (11.31) with critical exponent (II. 4p . One may regard the bound state iph as an orbit of (11.31) . 
From [TT], the orbital stability of iph is defined as follows: For all e > 0, there exists 6 > such 
that if ||'?/'o — Wftlli^i < 6 and -i/^ is a solution of (II. 3p in some interval [0, to) with ip\t=o = 'ipo, then 
ip{-, t) can be extended to a solution in < t < oo and supo<t<oo inf^eR ||'^('5 t)—4'h{-, s)\\h^ < 
Otherwise, the orbit iph is called orbital unstable. 

The functions V = V{x) and m = m(x) may play a crucial role on the orbital stability of 
iph- When m = 1 and V is of class (y)a and fulfills other conditions in [2n]-[3D], the orbital 
stability and instability of iph for the equation (II. 3p was established by Lin and Wei [23] if V 
has non-degenerate critical points. Under different conditions, e.g., h = 1 and A is large, results 
of the orbital stability problem can be found in [15|. One may also remark that the orbital 
stability problem of NLS with inhomogeneous nonlinearity has been investigated in [5] but only 
for the subcritical case, i.e., 1 < p < 1 + 

To state our main results, we need to introduce some notations. It is well-known that the 
positive solution of 

w + w^ = in , 

= max w(y) , wiy) — )■ as \y\ — +oo . i^-^) 

is radial [16] and unique [21]. We denote the solution and its linearized operator a.s w = w{r) 
and 

Lo := A - 1 + pwP-\ (1.7) 
respectively. For the orbital stability of tph, we set 

Lh:=h^A-{V + X) + mpul-^ (1.8) 

as the linearized operator of (II. 5p with respect to Uh and 

'H^ .o 1 o 1 




d{\) 



dx, (1.9) 



as the energy of Uh- Observe that Uh may depend on the variable A. Assume that d{X) is 
non-degenerate, i.e., d"{\) ^ 0. Let p{d") = 1 if rf" > 0; p{d") = if rf" < 0, and n{Lh) be the 
number of positive eigenvalues of L^. According to general theory of orbital stability of bound 
states (cf. [T7], [IS]), iph is orbital stable if n{Lh) = p{d"), and orbital unstable if n{Lh) —p{d") 
is odd (see page 309 of [IB])- It is remarkable that if both V and m are constant and p = 1 + -^, 
then d"{X) = 0. Consequently, from now on, we consider the critical exponent p = 1 + and 
assume the point Xq as a non-degenerate critical point of the function G defined by (cf. [20] , 
[39]) 

G{x) := [Vix) + A] m-^/\x) , Vx G , (1.10) 

provided V ^ and m > in M^. When ^ = in M^, xq is set as a non-degenerate critical 
point of the function m. 

For simplicity, we firstly switch off the potential V and obtain the following result. 
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Theorem 1.1. Let N < 3 be a positive integer, p = 1 + and the potential V = 0. Assume 
the function m = m{x) satisfies 

m G C^(M^);0 < tuq < m{x) < mi < oo; \m^'\x)\ < Cexp(7|x|), i = 1,2,3,4, (1.11) 

where mo, mi, 7 and C are positive constants, and m^'^\x) are the i-th derivatives of m{x). 
Suppose also that xq be a non- degenerate critical point of m[x) (xq is independent of X). Let 
il}h{x,t) := e'^^^'^Uhix) be a bound state of ( fi.51) . where Uh is a single-spike solution of ( fi.5|) 
concentrating at xq. Assume also 

m{xo)A^m{xo) < CN,i\Am{xo)\'^ + Cn,2 N\\V^m{xo)\\l - \Am{xo)(' 

+Cjv,3m(xo)V(Am)(xo) ■ [V'm(xo)] ^ V(Am)(,xo) , (1.12) 

where 



00 



2{N + 2)2 / v^+^wPLq^ {r^wP)dr 
Cn,i = '-^ , (1.13) 





00 



A{N + 2) J r^+^wP%dr 



Cn,2 = ^ , (1.14) 

iV2 J r^+^njP+^dr 


(iV + 2) ( / r^+i^f+Mr) 

Civ,3 = — 55 5S ) (1-15) 

jY j r^-iy^P+idr J r^+^vjP+'^dr 


are constants depending only on N. Here $0 = *l'o(^) satisfies 

^0 + ^^^^^0 - % + pwP^^<^o - ^$o-rV = 0,r = \x\e (0,oo), 
$o(0) = $'0(0) = 0. 

where Lq is defined in ( [i. 7| j. Then for any X > 0, iph is orbitally stable if h is sufficiently 
small and Xq is a non-degenerate local maximum point of the function m. Furthermore, for any 
X > 0, i/jh is orbitally unstable if h is sufficiently small and the number of positive eigenvalues 
of the Hessian matrix V'^m{xo) is odd. 

Remark 1: When = 1, xq = and the function m satisfies m"'(xo) = 0, (see (C.2) 
of [10]), the condition (11.121) of Theorem 1 1.1 1 is exactly same as the condition (4.14) of [10]. For 
N > 2, G.Fibich and X. -P. Wang (cf. [12]) considered the function m with radial symmetry, i.e., 
m = m{r),r = \x\ and m"'(0) = 0, and studied the orbital stability problem only for radial 
perturbations. Here we may include the case that the function m is not radially symmetric 
and the third order derivatives of the function m at xq can be nonzero. Moreover, we study 
the orbital stability problem for general perturbations including the non-radial perturbations. 
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Consequently, Theorem II .11 can be regarded as the most general theorem on the orbital stability 
problem of semiclassical NLS equations with critical exponent and nonlinear OLs. 

When the potential V is turned on, we may generalize the argument of Theorem 11.11 to 
obtain three theorems as follows: 

Theorem 1.2. Let N < 3 be a positive integer, p = 1 + Assume both the potential 
V = V{x) and the function m = m{x) satisfy the following conditions: there exist positive 
constants Vq, Vi, mo, mi, 7 and C such that 

V e C\R^y,0 <Vo <V{x) <Vi < 00; < Cexp(7|x|), z = 1,2, (1.17) 

and 

m G C^(M^);0 < mo < m(x) < mi < 00; |m(^)(x)| < Cexp(7|x|), z = 1, 2, (1.18) 

where V^'^\x) , m^'^\x) are the i-th derivatives of V (x) , m{x) , respectively. Suppose also that xq 
be a non- degenerate critical point of the function G defined in 1^1. 1C\) for fixed A > (xq may 
depend on \). Let ^h{x,t) := e*'^*/'^Uft,(x) be a bound state of lll.3\) . where Uh is a single-spike 
solution of ( (1.51) concentrating at xq . Then iph is orbitally unstable if h is sufficiently small and 
Xq is a non-degenerate local minimum point of G such that VV^(xo) 7^ 0. 

Theorem 1.3. Under the same hypotheses of Theorem M.^ assume also that W{xq) = and 
AV{xq) 7^ (thus Xq may be independent of X). Let n be the number of negative eigenvalues 
of the matrix V^G'(xo). Then iph is orbitally stable if h is sufficiently small and Xq is a non- 
degenerate local minimum point of G with AV{xq) > 0. Furthermore, iph is orbitally unstable 

if h is sufficiently small and n — | ^1 + j^^jf^ j is even. 

Theorem 1.4. Under the same hypotheses of Theorem assume also that 'WV{xq) = 0, 
AV{xq) = and U.ll\) holds for both V and m. Let n be the number of negative eigenvalues 
of the matrix V^G(xo). Suppose also that H{xq) > 0, where H{xq) defined in lji4-33\ ) involves 
the i-th derivatives (for < i < 4) of V and m at xq. Then iph is orbitally stable if h is 
sufficiently small and Xq is a non-degenerate local minimum point of G. Furthermore, iph is 
orbitally unstable if n is odd. 

Remark 2: Theorem ll.2til.4l may include all the cases of values VV^(xo) and AV{xq) for the 
orbital stability problem of (11.31) with critical exponent (II. 4p . Theorem 11.31 may generalize the 
main result of [23] to the case that the function m is a positive and nonconstant function. As 
= 0, Theorem II .41 coincides with Theorem II. II because of 

V'G(xo) = m(xo)~^^' [m(xo) V'y(xo) - y [V{xq) + A] V^m{xQ) 

Remark 3: In the following we give examples in dimension N = 2. Similar examples in 
dimension = 1 and 3 can also be given. Fist for x G M we define 

19 1 

Xi{x) = sin X H — sin^ x = - sin x sin(3x), 

6 8 24 

X2{x) = 2(1 — cosx) + -(1 — cosx)^ = o ~ o cosx + - cos(2x), 

3 2 3 6 

3 1 
X3(x) = sin'^x = - sinx — - sin(3x), 

X4(x) = 4(1 — cosx)^ = 6 — 8 cosx + 2 cos(2x). 
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respectively. Then Xi,X2,X3 and satisfy 

\Xi\ < l,X[{0) = l,xi'\0) = 0, for 2 < J < 4, 

< X2 < y,X^'(0) = 1,XJ')(0) = 0, for J = 1,3,4, 

iXsl < (0) = 1,X^,'\0) = 0, for J = 1,2,4, 

< X4 < 16, X[^\0) = 1, Xi^\0) = 0, for j = 1, 2, 3. 
Next for (x, y) G we set 

4 4 
V^(x, y) = ao + ^ aiXi(x) + + biX,{y), (1.19) 

i=l 1=1 

and 

4 4 

m(x, ?/) = Co + ^ QXi(x) + J2 d^X^{y), (1-20) 

i=l i=l 

where a^, fej, Cj, and di are constants. By the properties of Xi, X2, X3 and X4, the i-th derivatives 
of V and m at Xq = (0, 0) depend only on Oj, 6j and Cj, (ij respectively for 1 < i < 4. Recall that 
G{x, y) = [V{x, y) + A]m~^(x, y) for X = 2, we have 

VG'(O) = Co^(coai - (ao + A)ci,co6i - (oq + A)(ii)^, 



and then if VG'(O) = 0, 

V^G(O) = Co ^ 



coa2 - (flo + A)c2 

C062 - (flo + A)d2 



Now we can give examples for the potentials V and m which satisfy the assumptions in 
Theorems 1.2-1.4. 

(I) (Examples for Theorem 1.2) N = 2, Xq = (0, 0), V and m given in fll.l9p and fll.20p and 

ai,bi,Ci,di satisfy 

Co = flo + A, (fli, &i) = (ci, (ii) 7^ 0, 02 > C2 > 0, &2 > c?2 > 0, 

7 7 

and Co > 77(|ai| + |&i|),ao > Trdcil + |(ii|),ai = 6i = Ci = (ij = for z = 3,4, 
o o 

(II) (Examples for Theorem 1.3) First a special case for Theorem 1.3 is that Vm(xo) = 
0, V^m(xo) = and xo is a non-degenerate critical point of ^(x). Here we give another 
examples. The first one is in the stability case and the second is in the instability case. 

(a) (Stability) X = 2,xo = (0,0), V and m given in f ll.lQp and fll.20p and ai,bi,Ci,di 
satisfy 

32 32 

ao > 0, Co > — — C2 > 0, Co > — —d2 > 0, 02 > 0, 62 > 0, = 6i = Cj = rfj = for i = 1, 
o o 

then for any A > 0, the conditions in Theorem 1.3 for orbital stability will be 
satisfied. 
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(b) (Instability) N = 2,xo = (0, 0), V and m given in f ll.lQp and (ll.20p and Oj, 6j, q, di 
satisfy 

16 16 

ao > — —62 > 0, Co > — —C2 > 0, 02 + 62 > 0, d2 > 0, 
o o 

and ai = bi = Ci = di = for z = 1, 3, 4, 

then for any A > 0, the conditions in Theorem 1.3 for orbital instability will be 
satisfied. 

(Ill) (Examples for Theorem 1.4) Here we give two different examples. First we give examples 
in the case of 04 = 64 = 0. Specially, Theorem 1.4 is in this case. 

(a) (Stability) N = 2,xo = (0,0), V and m given in f ll.lQp and fll.20p and ai,bi,Ci,di 
satisfy 

32 32 

ao > 0,co > — — C2 > 0,co > ——d2 > 0, cq > — 32c4 > 0, |c2|, \d2\ small, or cq, |c4| large, 

o o 

and ai = hi = = = ds = d^ for i = 1,2, 3, 4, 

then for any A > 0, the conditions in Theorem 1.4 for orbital stability will be 
satisfied. Here \c2\, \d2\ small or cq, {c^l large are independent on A. 

(b) (Instability) N = 2,xo = (0, 0), V and m given in ( I1.19P and (ll.20p and a^, bi, Ci, di 
satisfy 

16 

ao > 0,Co > — — C2 > 0,(^2 > 0,co > — I6C4 > 0, |c2|, \d2\ small, or cq, {c^l large, 

and Qi = bi = = C3 = ds = d^ ioT i = 1, 2, 3, 4, 

then for any A > 0, the conditions in Theorem 1.4 for orbital instability will be 
satisfied. Here \c2\, \d2\ small or cq, |c4| large are independent on A. 

Second we give examples in the case of a4 + 64 7^ 0. 

1 (Stability) iV = 2,a;o = (0,0), V and m given in f irTOj) and ffL20|) and ai,k,Ci,di 
satisfy 

32 32 
ao > 0, Co > — — C2 > 0, Co > — —d2 > 0,a^> 0, 64 > 0, (a4 + 64) large, 

o o 

and aj = 6j = = C3 = C4 = = (i4 for 2 = 1, 2, 3, 

then for fixed A > 0, the conditions in Theorem 1.4 for orbital stability will be satisfied. 
Here (a4 + 64) large may depend on A. 

2 (Instability) = 2, xo = (0,0), V and m given in fll.lQp and fll.20p and ai,bi,Ci,di 
satisfy 

16 

ao > 0, Co > C2 > 0,d2 > 0, a4 > 0, 64 > 0, (a4 + 64) large, 

3 

and tti = bi = = cs = = ds = d^ for i = 1, 2, 3, 4, 

then for fixed A > 0, the conditions in Theorem 1.4 for orbital instability will be 
satisfied. Here (a4 + 64) large may depend on A. 



7 



The rest of this paper is organized as follows: In Section 2, we show the properties of Uh- 
Then we state the proof of Theorem 11.11 in Section 3. Theorem 1 1 . 2m .41 are proved in Section 4. 
Acknowledgments: The research of the first author is partially supported by a grant from 
NCTS and NSC of Taiwan. The research of the second author is partially supported by an 
Earmarked Grant from RGC of Hong Kong. 



2 Preliminaries 



In this section, we study the properties of a single-spike bound state of ( II. 5p concentrated 
at a non-degenerate critical point oiG{x) := \y{x) + A] m"^/^(x) (cf. [20], [39]). Let Xh be the 
unique local maximum point oi Uh- So — xq as /i — t- 0. 

Let Vh{y) '■= Uh{hy + Xh) for all y G M^. Then by (11.51) . Vh is a positive solution of 



At; - \V{hy + Xh) + X]v + m{hy + Xh)v^ = 0. 
For notation convenience, we still denote 

Lh:= A - [V{hy + xt) + A] + m{hy + Xh)pvl~^ 



(2.1) 



(2.2) 



as the linearized operator of the equation ( 12. ip with respect to the solution Vh- As the result 
of [39], Vh can be written as Vh = + where Wxf^ is the unique positive solution of 



Aw-[V{xh) + X\w + m{xh)w'P = Q in , 
ty(0) = max w{y) , w{y) — as \y\ — )■ +oo , 



and 

||0?i||oo as /i . 

Moreover, 

Vh{y) < C|2/|'^exp( - F'/'lyl) , Vy G , 
where V := inf^iv [V^(a;) + A] . From (12. 3p . it is easy to check that 

Wxu{y) = [V{xh) + X\^m{xh)~^w{^yV{xh) + Xy) 



(2.3) 

(2.4) 
(2.5) 

(2.6) 



where w is the positive solution of (II. 6p . 

For the single-spike solution of (II. 5p . we recall the following result from [38] and [39] : 

Lemma 2.1. Assume that there are positive constants 7 and C such that 

\VV{x)\, |Vm(x)| < Cexp(7|x|) , Vx G . 

Then 

-Vm{hy + Xh)vl^^ - - W{hy + Xh)vl dy = 



p + 



(2.7) 
(2.8) 



for < h < Hq, where Hq is a positive constant depending on 7 and A. 
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In the rest of this section, for simphcity, we switch off the potential V, i.e., set ^ = 0. Then 
by Lemma 12. 1^ we obtain the uniqueness of Uh as follows: 

Lemma 2.2. Suppose \2. 7\ ) holds, ^ = and xq is a non- degenerate critical point of m. Then 
Uh is unique. 

Proof. Suppose u^ i and m/j 2 are different single-spike solutions of (II. 5p concentrating at the 
same point Xq. Let vi{y) := Uh,i{hy + Xq) and f2(y) := Uh^2{hy + Xq). Then both Vi and V2 
satisfy 

Av - \v + m{hy + Xq)v'p = , for y G , 
and f 1, f 2 — uniformly on as /i — t- 0. Due to Vi ^ ^2, we may set 

Vi - V2 
Vh ■= J, 77 , 

and then Vh satisfies 

Avh - Xvh + m{xo)pw^~^Vh + [Tn{hy + a;o) - m{xo)]pwP~^Vh + N{vh) = 0, (2.9) 

where N(vh) = rn{hy + xo)[vi—V2—pw^~^{vi — V2)]/\\vi — V2\\oo- Hence by the standard elliptic 
PDE theorems on the equation 02. 9p . we may take a subsequence Vh -^Vo, where Vq solves 

Avq -vo + m{xo)pw^~^vo = 0. 

Consequently, there exist constants c/s such that 

N 

vo = ^ CjdjWxQ . (2.10) 
i=i 

Let yh be such that Vh{yh) = \\vh\\oo = 1 (the same proof applies if Vh{yh) = — !)• Then by 
the Maximum Principle, we have \yh\ < C. On the other hand, as (12.81) . we may obtain 

Vm{hy + Xo)f i^^rfy = = j Vm{hy + XQ)v2^^dy. 
Thus 

/p+i p+i , 
Vm{hy + xo) n '""^ )vhdy = 0. (2.11) 
V V1-V2 / 

Note that for alH = 1, ■ ■ ■ , A^, as /i — t- 0, 

dim{hy + xq) = hS^ dikm{xo)yk + o{h) , and — = (p + l)w^„ + o(l) . 

Hence from (12.101) and (12.111) . we may obtain 

N N 



■"^xo jdy + o{h) 

jjJV k=l j = l 

N 

= -h^dijm{xo)cj / w^^^dy + o{h). 

i=i j.r 



Hence by the assumption that V^m(xo) is non-degenerate, Cj = for j = 1, ■ ■ ■ , A^, i.e., vq = 0. 
This may contradict to the fact that 1 = Vh{yh) ~^ '^o(yo) for some yo G M^. Therefore, we may 
complete the proof of Lemma 12.21 □ 



By Lemma 12.11 we may simphfy the proof of [21] and get a shorter proof of the asymptotic 
behavior of Xji S clS follows: 



(2.12) 



Lemma 2.3. Under the same hypotheses of Lemma\K 



Xh = + o{h) as h . 



Proof. Fix i G {1, ■ ■ ■ ,N} arbitrarily. By Taylor's expansion of dim{x) and Vm(xo) = 0, we 
obtain 



N 



dim{hy + Xh) = ^dijm{xo){hyj + Xhj - xqj) + o{h) + o{\xh - xo\). 
Hence by Lemma [2.11 and Vh = + o(l), we have 







dim{hy + Xhjv^j^^dy 



N 

^dijm{xo){xhj - xoj) / w^-^^dy 

i=i Ar 



+ o{h) +o{\xh-xo\) 



Here we have used the fact that J yjW^'^^ dy = for j = 1, ■ ■ ■ , N. Using the assumption that 
V^m(xo) is non-degenerate, we obtain f l2.12p . □ 



Following the idea of [25], we may use Lemma [2.31 to show the asymptotic behavior of as 
follows: 



Lemma 2.4. Under the same hypotheses of Lemma \2.^ 

Vh = uix^ + /i^02 + o(/i^) , as h ^ , 



(2.13) 



where 02 satisfies 

1 ^ 

A(/.2 - X4>2 + m{xh)pwP;^(P2 9-M''^)yiy'j<n = , and V</.2(0) = 0. (2.14) 

Proof. Let (ph = Vh — Wx^- Then it is easy to check that |0ft,| — )■ uniformly, and (ph satisfies 
- + m{hy + x;,)p<- V/. + iV(0h) + R{<Ph) = , and V0;,(O) = 0, (2.15) 

where 



and 



N{(t)h) = m{hy + Xh 

R{4)h) = m{hy + Xh) - m{xh) 



10 



Note that by Lemma [2.31 and Vm(xo) = 0, 



m{hy + Xh) - m{xh) =hy ■ Vm{xh) + y X] (^ij^i^hiyiVj + o{h^ 

= — ^ dijm{xQ)yiyj + o{h^). 



Now we claim that \(f)h\ < ch"^ by contradiction. Suppose that h'^ 



-'h oo 



(2.16) 



— !■ oo. Let 



/\\<Ph\\oD- Then (j)h satisfies 



(2.17) 



Note that by (HTT^ . 



||</'h||c50 II </'/i||oo 



(2.18) 



Let yh be such that (phiVh) = ||0/i||oo = 1 (the same proof apphes if iphiVh) = — !)• Then by 
(12.171) — (12. ISp and the Maximum Principle, we have \yh\ < C. On the other hand, by the usual 
elliptic regularity theory, we may take a subsequence (ph ^ (po, where 0o satisfies 

A0O - A0O + m(a;o)p<;Vo = , and V0o(O) = 0. 

Hence 0o = 0. This may contradict to the fact that 1 = (phiVh) 4>oiyo) for some yo. Therefore, 
we may complete the claim that \(j)h\ < ch"^. 

Now we set (l)h,2 = 4>h — h'^4>2- Then (f)h,2 = 0(/i^) and satisfies 



Mh,2 - Mh,2 + m{hy + Xh)pwl-^<Ph,2 + N{(Ph^2) + R{<Ph,2) = , and V(/.,,,2(0) = 



where 



and 



N{4>h,2) = m{hy + Xh) {w,^ + K'ct)^ + <i)h,2Y - - pwP;\h^(P2 + 4>h,2) 



R{4>h,2) = m{hy + Xh) - m{xh) dijm{xo)yiyj 



m{hy + Xh) - m{xh) 



Thus as for previous argument, we may have (j)h,2 = o(/i^) and complete the proof of Lemma [231 



□ 



As for Proposition 3.1 of [23], one may get two lemmas as follows: 
Lemma 2.5. For h small enough, the maps 

L^J := A0 - [V{xh) + A]0 + m(a;/,)p<^V 



11 



are uniformly invertible from K^^^ to Cj; , where 



K. 



(0Gi/2(M^) f ^d,w^^Jy = Q,3 = l,--- ,n\ dH 



(pdjWr^^dy = 0, j 



Lemma 2.6. The map 



has eigenvalues fij , j 



- A0 - [V{xo) + A] + m{xo)pwP^ V 
■ ■ ■ , + 2 satisfying 

fll > = fl2 = ■ ■ ■ = fJ'N+l > fJ'N+2 > ■ ■ ■ , 

where the kernel of L^^ is spanned by djWx^, j = 1, ■ ■ ■ , N and fii is simple. 

In this section, our main result is the small eigenvalue estimates of Lh given by 



Theorem 2.7. Under the same hypotheses of Lemma \2.Sl for h small enough, the eigenvalue 
problem 



has exactly N eigenvalues ix'f^ , j 



,N, in the interval which satisfy 



[A 



CqU. 



[up to a subsequence) as h ^ , for j = 1, 



AT, 



(2.19) 



(2.20) 



where Hi and ij.n+2 o'^e defined in Lemma [27^ uj's are the eigenvalues of the Hessian matrix 
V^m(xo) and Cq = 2m{xo) ^ positive constant. Furthermore, the corresponding eigenf unctions 
ipj^ 's satisfy 



N 



where Uj = (a^ 



, aNj) 



oCDla-m 



3 



AT, 



(2.21) 



1=1 



is the eigenvector associated with Uj, namely, 

V^m(xo)aj = ujaj . (2.22) 

Here o(l) is a small quantity tending to zero and 0(1) is a bounded quantity as h goes to zero. 

Remark 4: (1) Since converges to L^^ in the strong resolvent sense, in the interval (|/ii, oo) 
Lh has only one positive eigenvalues which is simple and goes to jii as h goes to 0. 

(2) After changing variables t t/h,y = {x—XQ)/h, L^ becomes —Rh, which is the notation 
used in page 190 of [17]. Thus the number of negative eigenvalues of Rh equals the number of 
positive eigenvalues of Lh, which we denote by n{Lh). 

(3) By (12.201) . the sign of small eigenvalue /x;^ of Lh is the same as the one of eigenvalue 
Uj of V^m(xo). If we denote the number of positive eigenvalues of V^m(xo) by n, then the 
number of positive eigenvalues of Lh in the interval \ fJ'N+2] equals n. Adding another one 
in the interval (|/ii, 00), the number of positive eigenvalues of Lh equals n + 1. In particular, 
if V^m(a;o) is negative definite, then n = and thus n{Lh) = 1. 
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Proof. We may follow the arguments given in Section 5 of ^2]. Assume that ||(y9/i||/^2 = 1. By 
Lemma (12.61) it is easy to see that /x/i — )■ as /i — )■ 0, where fih ^ {/^L ■ ■ ■ ^l^h}- Then the 
corresponding eigenfunctions iphS can be written as 



N 



where (/?^ G K^^. Hence by (ICTD and fl2:23|) . </?^ satisfies 



(2.23) 



(2.24) 



where 



) = "^{hy + Xh,)p(i;^ ^ - wl^ + m(/i?/ + Xh) - m{xh) 

Using (12.161) and Lemma [2.41 we have 

LhdjW^^ = m{hy + Xh)p{vl"^ - wl~^)djW^^ + m{hy + x/^) - m{xh) pwl'^djW^^ = 0{h^). 

(2.25) 

From Lemma [2751 the map L^^ = A — \ + m{xh)p'w^~^ is uniformly invertible in the space K^^. 
Thus by (I2.25P and ///i — j- 0, we have 



N 



yi\\H^<c{h^ + \fih\)J2 



K\ ■ 



(2.26) 



To estimate fih and a-^'s, multiplying (I2.24p by dkiv^,^ and integrating over , we may 
obtain 

N N 

{Lhipi) dkW^,dy + y2^h {LhdjW^J dkW^^dy = /i/i / djW^,^dkW^,^dy . (2.27) 

Here we have used the fact that G -ft'j"^. Using (12.251) . (12.261) . = o(l) and integration by 
parts, we obtain 



and 



h 



{LhdjWxJ dkW^^dy = — — - / w^'j^^dydjkm{xo) + o(/i^) 



(2.28) 



(2.29) 



P + 1 jRiV 

which we have proved in Appendix A. Substituting (12.281) and (12.291) into (12.271) . we may obtain 



P 



1 f ^ 

— T / Kt^dyy^djkmixo) 



, _ l^h k 



I / [dkW^.Ydy + o{l). 
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Since ||v5fc||L2 = 1, f !2.23p implies that a/j := (a^, ■ ■ ■ ,a^)'^ is bound. Moreover, by fl2.26p . 
does not converge to 0. Thus p- — )■ CoVj ioi j = 1, ■ ■ ■ , N and ah — )■ aj , where 



N 



Co 



[p+l) J^N \Vwcco\^dy 2m(xo)' 
and aj is the eigenvector corresponding to uj. Here we have use the fact that 



/ \Vwa:o\'^dy 



A^ + 2 



m(xo) / vjPj dy 



which can be proved by Pohozeve identity. The rest of the proof follows from a perturbation 
result, similar to page 1473-1474 of j42]. We may omit the details here. □ 



3 Proof of Theorem 11.1 



In this Section, we firstly study the asymptotic expansion of d"{X) as /i — j- 0, and then 
complete the proof of Theorem 1 1.1[ To drive the 0{h^) order terms of d"{\)/h^ , we need the 
following lemma: 



Lemma 3.1. Under the same hypotheses of Lemma \2.Br 

Xh = Xq + h'^Xi + 0{h^) , as /i — )■ , 

where Xi G satisfies 

J lyl'^w^'^^dy 

Proof. By Lemma 12.31 and Vm(xo) = 0, for alH = 1, ■ ■ ■ ,N, we have 

N 

dim{hy + x^) = ^ dijm{xo) {hyj + Xh,j - XQ,j) + 0{h^ 
i=i 

Then by (12.81) , (13. 3p and Lemma 12. 4[ we have 

p+i 



= / dim{hy + Xh)vl dy 



N 



dijm{xo) J {hyj + Xhj - xqj) + 0{h) dy + 0{h^ 

N 

^ dijm{xo) {xh,j - Xqj) / w^^^dy + 0(/i^) . 



(3.1) 



(3.2) 



(3.3) 
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Here we have used the fact that / UjW'^:^^ dy = for j = 1, ■ ■ ■ , A^. Thus Xh = Xq + 0(/i^). 

Consequently, we may set Xh = Xq + h?Xh- Then Xh = 0(1) and by Taylor's formula of dim{x), 
we have 

dim{hy + x?,) = ^ dijm{xo) {hyj + h'^Xh) + y X] ^ufc"^(^o)z/jl/fc + 0{h^). (3.4) 
j=i j,k=i 

Hence byf l2.8l) . f l3.4p and Lemma [2.41 we may obtain 

N „ , 9 N 



=h'^^dijm{xo)xh,j J wl'^^dy + y ^ 9ijkm{xo) J yjykwl'j^^dy + 0{h^ 
^ f ^ f 



Here we have used the fact that 

Jy,wP^'dy = 0, Vj = l,---,iV, 

/ = % / , Vj, A: = 1, • ■ ■ , AT. 

Therefore, we may complete the proof because 

w^xo(2/) = A^/^m(xo)"^/^^(yA2/). 



□ 



From Lemma 12.41 and 13.11 we may deduce that 

Theorem 3.2. Under the same hypotheses of Lemma \2.2\ for h small enough, Uh is smooth 
on \. Let Rh := ^{hy + Xh). Then 

LhRh - Vh = 0. (3.5) 

and 

N 

Rh = Ro + Y, 4djW,, + h^Ri + Ri , (3.6) 
i=i 

where Rq = X^^{^;ziVh + \y ■ Vi^h), = 0{h), Rj[ = 0{h^) and Ri satisfies 



ARi - XRi + m{xh)pwl^ "^^i ~ ^ X] ^u"^(^o)yi% w^Xh = ■ (3-7) 

Furthermore, 

f \y\'^w^^^dy 

V'm{xo){h~'c,) ^ -^^^j-^^ViAm){xo) , as h^O, (3.8) 

R^ 

where := (4, ■ ■ ■ ,c^)'^. 
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Proof. By Lemma [2.21 and Theorem 12. 7[ Uh is unique and non-degenerate. Consequently, Uh is 
smooth on A and satisfies fl3.5p . Now we decompose Rh as 



N 



h 5 



where Rj^ G K^^. Then Rf^ satisfies 



N 



LhRi + [LhRo + h^LhRi -vyi\+Y^ ^LhdjW^^ = 0. 



As for the proof of Theorem 12.71 we have 



N 



It is easy to check 



\Rt\\H^ < c( \\LhRo + h'^LhRi - Vhh^ + J2 \^h\h' 



L/,i?o = Vh - ■ VmQiy + Xh)vl. 



Hence by Lemma [231 13. 1[ (13. 7p and ( 13. lip , we obtain 



LhRq + h LhRi — Vh 



2J = 1 

Consequently, by fl3.10p . 



i,j,k=l 



N 

\Ri\\m<c(^k' + J2\<\^' 



i=i 



To estimate c^'s, we may multiply (13. 9p by dkWx^ and integrate over M . Then 



{LhRi)dkWx^dy + 



-^/i-Ro + h LhRi — Vh 



dkW^^dy 



^ ■ f 



Hence by ([2:291) . iKWt may imply 
C 



K\< 



{LhRi)dkWx^dy\ + I / [LhRo + h^LhRi - Vh\dkWx^dy\ 
Jrn 

Using integration by parts and (I2.25p . we have 

/ {LhRi)dkWx^dy = j Rj^LhdkW^^dy = \\Rj;\\L20{h?') . 



(3.9) 



(3.10) 



(3.11) 



~ 2A [ ^ 9ijm(a;o)xi,iyj + 2 5Z 9ijkm{xo)yiyjyk w^^ + 0(/i^) . (3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 
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Therefore, by fl3.12p . fl3.13p . (13.151) and (13.161) . we may obtain |c^| = 0{h). Consequently, 
by dSm, Ri = 0{h^). Thus by dsn. 



{LhRi)duw^Jy = 0(/i^). 



(3.17) 



Hence by dOHl), dffl) and fIXTTD . flXTI|) gives 



1 



N 



~^dy^djkm{xo)[h ^4) 



TV 



2A /^». I X] ^^j"^*^^o)^i.»%' + 9 dijim{xo)yiyjyi wl^dkW ^^dy + o{l). (3.18) 



Using integration by parts, we obtain 



JgAT 

/ yiVjViwP^dkW^Jy : 



N{p+1) 



I \y?<'dy, 



where 6 is the Kronecker symbol. Hence by (I3.18p . |c^| = 0{h) for j = 1, ■ ■ ■ , N. Moreover, 
by (13. 2p . we obtain (13. 8 p and complete the proof. 

Let us now compute d"{\). From (II. 9p . it is easy to get 



d'{X) 



u\dx 



□ 



and hence 



d"{\) = J Uh^^dx = J ^hRhdy ■ 



(3.19) 



Using integration by parts and (13. 5p . we have 

//"II 1 N r 

VhRody = j Vh\~\j—^Vh + -y ■ Vvh)dy = A'^^ - J ^Idy = , (3.20) 



since p = 1 + Hence, by (I3.19P and Theorem 13.21 we have 



N 



Vh 



Ro + Y^ 4p3^x^ + h^Ri + K dy 



N 



Vh 1^ cldjWxf^ + h?Ri + Rh dy (because j VhRody = O) 

N 

Rh [ ^ (^h^hdjWxh + h'^LhRi + ivfti?^ (because L/,i?/i = 

Af AT 
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Therefore, by 1^2^ . ( K9\\ and = 0{h). 
d"{X) 



N 



Rn 



Vh — LhRo 



j,k=i 



dy+ / dkW^^{LhdjW^^)dy 



+ h' j Ri{URi)dy + 0{h'). 



(3.21) 



For the integral / Rq 



dy, by (13. Ill) and using integration by parts, we have 



Rn 



Vh — LhRo 



/I 1 r /i 

A"^ {^^Zrj^h + ■ ^^Vh) ^y ■ Vm{hy + Xh)vl dy 



N 



N 



2A2 J 4(iV + 2) 
Note that by Lemma 12.41 13.11 and Theorem 13. 2[ we have 

N 

hy ■ VmQiy + Xh) - X] ^v^i^V + ^h)yiyj 



hy ■ VmQiy + Xh) - X] dijUiQiy + Xh)yiyj 



vl^^dy. 



--hy ■ Vm{xh) - y X] ^ijkm{xh)yiyjyk - — ^ dijkim{xh)yiyjykyi + o{h'^) 



N 



2 ^ ' ' '^^'^ 3 

i,j,k=l i,j,k,l=l 



and 



Hence 



VI. = w 



(3.22) 



Rq 



Vh — LhRo 



dy 



N 



3 X] 9ijkim{xh)yiyjykyi w^+^dy + o(/i^) 

i,j,k,l=l 

^/^^^. A-^xo)"^'^ I lyl'w^-^'dyA'mixo) + o{h') . (3.23) 



Here we have used the following identities: 

/ yiW^t^dy = J yiyjykwl+^dy = , for all i, j, = 1, ■ ■ ■ , ; 

/ yiyjykyi'^xt^dy = O , if yiyjykyi is an odd function on one of its variate ; 

7 yf<'dy = jrmi) I \y\'<'dy , for aiH = i, ■ ■ ■ , tv ; 

RN 

I yh]<t^dy = nWW) I \y\^<t^dy , for all i j , 
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which can be proved by polar coordinates. 

N 

For the sum ^ c;^c| J dkWx^{LhdjWx,^)dy, we may use fl2.29p and (13. 8p to get 

j,k=l MAT 



^ • f 



j,k=l 



P+l 
h 



j,k=i 



( / \y\'w^+'dyy 



(3.24) 

For the integral J Ri{LhRi)dy, by (13. 7p . it is obvious that Ri{\^^y) satisfies 

1 iV iv ^ 

AR-R + pw^-^R - -A^-2m(x/,)^T-i ^ dijm{xo)yiyjwP = . (3.25) 

Hence 

/i^y R^{L,,R,)dy = j i?i(L,.,i?i)rfy + 0(/i6) 

=— A"'^m(xo)"^"^ ^ dijm{xo)dkim{xo) / yiyiw'^L'^^{ykyiw^)dy + 0{h^) 
A"=^m(xo)~^~'|Am(xo)|' /" r^w^ L^\r''w^)dy 



4Ar2 

+ ^^^^^^^ A-^m(xo)"^"lV^m(xo)||^ j r^w^%{r)dy 



4 



2A^2(]v + 2) 



X-^m{xo)-^-^\Am{xo)\^ j r V<l>o(r)rf?/ + 0(/i*^). (3.26) 
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N 

Here || V^m(a;o)||2 = Yl ''^Iji^o) and we have used the following identities: 

ylw^L^\yW)dy = ^ | r'w^L^\r'wndy + ^^^^^^ / r'w^Mr)dy , (3.27) 

RJV 



yN^WNW^L^ ^{yN^iyNwP)dy = ^ / r^u7*'$o(r)rfy , (3.29) 



where $o satisfies (11.161) . which we have proved in Appendix B. 
Therefore, combining ([MID, <^M>, (13:21) and (KW> . we obtain 



A-3m(xo)-^~' y li/l V+M?/A'm(xo) 



8(A^ + 2)2 



+ 8iV(^\2) ^"'"^^^°^""" V(Am)(xo) ■ [V^m(xo)] " V(Am)(a:o) 

R^ 

+ ^X-'m{xo)-^-'\Am{xo)\' I lyl'w^L^' {\y\'wP)dy 

R^ 

+ ^^^^^^^^^^A-=^m(xo)-^-'[iV||V2m(xo)||^-|Am(xo)pJ J \y\^wP%{\y\)dy. 

R^ 

Consequently, 

8(iV + 2)^m(a;o)f+2A 

RJV 



-ci"(A) =C;v,i|Am(xo)P + C^,2 (iV|| V2m(xo) ||^ - |Am(xo)P) 

+ CN,-Mxo)[v{Am){xo) ■ [V'm(xo)]"V(Am)(xo) 

— m(xo)A^m(xo) + o(l) , 

where Cn,i,Cn,2,Cn,3 are constants given by f ll.l3p . (11.141) . (ll.lSp . respectively. 

Now we may prove Theorem 11.11 as follows: Suppose that xq is a non-degenerate local 
maximum point of the function m{x), then the Hessian matrix V^m(xo) of m at xq is negative 
definite. By Theorem 12.71 we have n{Lh) = 1. On the other hand, we have p{d") = 1. Thus iph 
is orbital stable by the orbital stability criteria of [I7]-[l8]. For orbital instability, we denote 
the number of positive eigenvalues of the Hessian matrix V^m(xo) by n. Then by Theorem 12.71 
we obtain n{Lh) = n + 1. On the other hand, we have p{d") = 1. Thus by the instability 
criteria of [18], we conclude that iph is orbitally unstable if n is odd. This may complete the 
proof of Theorem 11.11 



20 



4 Proof of Theorem ll.2til.4 



In this section, we may generalize the argument of Section 2 and 3 to prove Theorem ll.2til.4[ 
Let Vh{y) ■= Uhihy + Xh), where Uh is a single-spike bound state of (II. 5p with a unique local 
maximum point at Xh. Then Vh satisfies 



V{hy + Xh) + \ Vh + m{hy + Xh)vl = in 



Suppose (ETD hold. By (ESD and [39], we have 



N 

m(xo)Vy(xo) = — [V{xq) + A] Vm(xo) 



(4.1) 



(4.2) 



so Xq may depend on A. Note that by (14.21) . Vm(xo) = if and only if VV^(xo) = 0. By direct 
computation on the function G, 

SyG(xo) = m(xo)"^"^ m{xQ)dijV{xQ) + (1 - ^)diV{xQ)djm{xQ) 

N 



In particular, if Vm(a;o) = 0, then 
V^G'(xo) = m(a:o)"^^ 



mi 



-[V(xo) + A]9.,m(xo) 



;a:o)VV(xo) - y [^^(xo) + \\V^m{xo] 



Using the identity (12. 8p . one may follow the arguments of Lemma \2 . 2112 .41 to get the unique- 
ness of Uh and 



Xh =XQ + o{h) ; 

Vh = + h(t)i + /j^02 + o(/i^) , 
where 0i and 02 satisfy V</)i(0) = V02(O) = , 

A01 - [V"(xo) + A] 01 + m(xo)ptfS,7 Vi - y ■ W{xq)wxo + y ■ Vm(xo)tf^Q = 

and 

N 



(4.3) 
(4.4) 



(4.5) 



A02 - [V"(xft,) + A]02 + m{xh)pwl^ V2 - y ■ VV"(xo)0i " ^ X] dijV{xo)yiyjW^^ 

1 ^ 1 

+ y ■ Vm{xo)pwl^ ^'^i + 9 X] dijm{xo)yiyjwP^ + -m{xo)p{p - l)w^~^(pl = . 



(4.6) 



Here we have used the hypothesis that xq is a non-degenerate point of the function G. And the 
only difference in the proof is that we need to estimate the term 



P 



^Vm(xo) j vl'^^dy -^W{xo) j vldy , 
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to estimate which one may use the following Pohozaev identity (cf. [M] ) 



h 



N + 



-m{hy + Xh) H —y ■ Wm{hy + Xh) v'h dy 



h 



V{hy + Xh) + \ + -y W{hy + Xh) v^dy . 



For the small eigenvalue estimates of L^, one may generalize the idea of Theorem 12 .71 to get 
Theorem 4.1. Fork small enough, the eigenvalue problem 

Lh<^h = fJ'h'^h (4.7) 



has exactly N eigenvalues fi^l ,j = 1, - ■ ■ N, in the interval |/iAf+2]; which satisfy and 



-j^^coUj, as h^O, forj = l,---N, 



(4.8) 



where fii and fiN+2 o,i"^ defined Lemma \2.6l Uj's are the eigenvalues of the Hessian matrix 
V^G(a;o), and cq = ~^y^jy^ = —G{xq)~^ is a negative constant. Furthermore, the corre- 
sponding eigenf unctions (pj^ 's satisfy 



N 



(4.9) 



i=l 



where each ipi is the solution of 



V{xh) + \ + m{xh)pwl~^'4)i 
y ■ W{xh) + y ■ Vm{xh)pwll^ + m{xh)p{p - l)ty^~Vi 
and aj = {aij, ■ ■ ■ , aNj)'^ is the eigenvector corresponding to uj, namely, 

V^G'(xo)aj = h'jaj . 



(4.10) 



(4.11) 



Remark 5: (1) To prove it, one may follow the arguments in the proof of Theorem 12.71 and 
use the following identity 



dkW^^^Lh{djW^^^ + hi)j)dy 



N + 2 



w^'+^dydjkGixo) + o{h' 



(4.12) 



to replace fl2.29p (see Appendix C). The main difference between Theorem 12.71 and 14.11 is that 
(14. 9 P has the solution tpi of f l4.10p which comes from 



LhdiWxt, = h -y ■ W{xo) + y ■ Vm{xo)pw^^ ^ 

+ m{xo)p{p - l)u;^;;Vi diWx^ + 0(/i^) 



(4.13) 
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(2) Let n be the number of negative eigenvalues of the matrix 5'^G{xq), then similar to the 
Remark 4(3), the number of positive eigenvalues of equals n + 1, i.e., n{Lh) = n + 1. 

Since the potential function V is nonzero, then xq may depend on A and the asymptotic 
expansion of d"{X) becomes more complicated. Indeed, when m = 1 and AV{xq) ^ 0, the 
result in [25] shows that the effect of potential function V on d"{X) is 0(/i^). On the other 
hand, when V = and condition f ll.l2p holds, the effect of m on d"{X) is 0(/i^) (see Section 3). 
Generally, when both m and V are not constant, we may show 

(I) The effect of V and m on d"{X) is 0(1) if VV{xo) ^ (see Theorem O]) ; 

(II) The effect of V and m on d"{X) is 0{h^) if VV{xo) = and AV{xo) ^ (see Theo- 
rem [L3]); 

(III) The effect of V and m on d"{X) is 0{h^) if VV{xo) = 0,AV{xo) = and some local 
condition hold (see Theorem ll.4p . 

Now we divide three cases to prove these results. 
Case I: VV{xo) ^ 0. 

Let Rh := ^{hy + Xh). Then ^ and ^Mj hold. Hence one may apply the idea of 
Theorem 13.21 to get 



N 




(4.14) 



i=l 



where as /i — > 0, c/i 



, ■ ■ ■ , c^) satisfies 



V^G{xo){hch) ^ 



m(xo) 2 ^Vm(xo) , 



(4.15) 



and 



Ro = [V{xh) + A]"'( -vh + -y ■ Vvh) , Ri = 0{h) 



(4.16) 



Thus 






N 



23 



Therefore, by flCTj) . ( 1413|) . (gH]), f lTO]) and dUe]), we obtain 
d"(A) 



4(iV + 2 



~N~2 



j wP-^^dy\/m{xo) ■ [V^G(xo)] Vm(xo) + 0(/i) . (4.17) 



Consequently, if xq is a non-degenerate local minimum point of G, then the Hessian matrix 
V'^G{xq) is positive definite. By Theorem 14.11 we have n{Lh) = 1. On the other hand, 
by f l4.17p . we have p{d") = 0. Thus we complete the proof of Theorem 11.21 by the orbital 
instability criteria of [T7]-|18j. 

Case II: VV{xo) = and AV{xo) ^ 0. 

Firstly, note that in this case, (pi = and ipi = 0. Then one may apply the idea of Lemma [XT] 
and Theorem 13.21 to obtain 

Xh =xo + + 0{h^) ; 



N 



Rh =Ro + ^ 49jW.H + h^Ri + Ri 



where xi G satisfies 



V^G(xo)xi = - [V{xo) + A] ^'mixor^^ 



f j \y\^w'^dy\ 



AN 



1 / N 



/ / \y\'^wP+^dy\ 



\ 



f wP+^dy 



f wP+^dy 
V(Am)(a;o), 



V(Ay)(xo) 



Ri satisfies 

Ai?i - [V{xh) + \\Ri + m{xh)pwl-^Ri 

N 

r , , -, _1 r ' 

ij=l ~ ij=l 

Rj[ = 0{h?) and = 0{h) for j = 1, ■ ■ ■ , N. Moreover, Ch := {c\, ■ ■ ■ , c^) satisfies 

V^Gixo){h-'c,) =co + o(l). 



[V"(xft,) + A] ^[Y^ dijV{xo)yiyjW^^ ~ 2 ^ 9ijm{xo)yiyjwP^ 







where 



Co = - [V{xo) + X] ^m(xo) ^VV(xo)xi 

/ / \y\Vdy\ 



N + 2 
2N 



[V{xo) + A] m{xQ) 2 



v 



/ wP+^dy 



V(Ay)(xo) 



\[V{xo) + \] ^m{xo)-^-^ 



I J \y\'^wP+^dy\ 

IN 

/ wP+^dy 



V{Am){xo). 



(4.18) 
(4.19) 



(4.20) 



(4.21) 



(4.22) 



(4.23) 
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Hence 



d"iX) 



VhRhdy = / Vh 



N 



Ro + Yl (^hdjW^n + h^Ri + Rh\ dy 
y w/j 1^ ^h^j'^^h + h^Ri + R-h dy ^because J v^R^dy = 
/ -Rft 1^ ^ c{LhdjW^^ + h^LhRi + LhRh dy (^because LhR, 



N 



N 



N 



Ro + Yl ^hdkWx^ + h^Ri + Rh\[Y^ ^LhdjW^^ + h'^LhRi + URt dy . 



k=l 



Therefore, by flCTj) . f Hl3|) and f lCTj) . we obtain 



d"i\) 



N 



Ro[vh - LhRo]dy + ^ 44 / dkWx^Lh{djWx^)dy 



j,k=i 



+ j Ri {LhRi)dy + 0{h^) . 



(4.24) 



Vh 



For the integral / Rq \vh — L^Rq^ dy, by direct computation, we have 

Vh - LhRo = - [Vixh) + A] ^Vihy + x^) - Vixh) + ^y ■ WV^hy + Xh) 
+ ^\^{^h) + y^] ■Vm{hy + Xh)vl. 
Thus by g2D, fl4A8|) and ([21]), we obtain 

R^[vh-LhR^]dy=^[V{x^) + \]-''m{x^)-'^ j \y\^w^dyAV{x^) + 0{h^) . 



(4.25) 



(4.26) 



N 



For the sum ^ 44 / dkW,,{LhdjW^,)dy, by dHO]), (il3D and 4 = 0{h) for j 

j,A:=l Riv 



1, ■ • ■ , iV, we have 



J2 44 / dkW^^{LhdjW^^)dy = 0(/i^) . 



(4.27) 



Combining Km . IKT7\\ and KMl . we obtain 
rf"(A) 



2N 



[V{xo) + A] 'm(xo)-^ / VdyAF(xo) + 0(/i^) . 



(4.28) 



25 



Consequently, by f l4.28p . we have p{d") = |(1 + On the other hand, by Theorem 14.11 

we have n{Lh) = n + 1. Thus we complete the proof of Theorem 11.31 by the orbital stability 
and instability criteria of [T7]-|18]. 

Case III: VV{xo) = , AV{xo) = 0. 

In this case, we shall use K23^ . flCTj) and (^2^ to compute the 0{h^) term oid"{X)/h^. 



For the integral J Rq 



Vh — LhRo 



dy, by (14.251) and integration by parts, we obtain 



Rn 



Vh — LhRo 



dy 



/IX h 
{j:ri"h + 2^ ■ ^^^^ [^^^^ + ~ ^^^^^ + 2^ ■ ^^^^y + ^h)\vhdy 

/II h 



8 



+ 



N 



[V{xh) + A] M ?,hy VV{hy + Xh) + h^Y^ di^V{hy + Xh)yiyj vldy 



N 



iVixH) + X]' 



N 



hy ■ Vm{hy + Xh) - dijinQiy + Xh)yiyj vl~^^dy . 



l{N + 2) 

Hence by (14. 180 . (I4.19P and Taylor's formulas of V and m, we have 



Rq 



Vh — LhRo 



dy 



-^[Vixh) + X\ 



N N 

Ah^ ^ dijV{xo)yiyjwl^ + Sh"^ ^ dijV {xo)yiyjW.^^(j)2 

i,j=l i,j=l 
N N 

+ Ah'^ ^ dijkV{xo)xi^iyjykwl^ + dijkiV{xo)yiyjykyiw, 

i,j,k=l i,j,k,l=l 



+ 



N 



8{N + 2) 



[V{xh) + X\ 



N 



^ dijkim{xo)yiyjykyi 



dy 



(4.29) 



i,j,k,l=l 



N 



For the sum ^ c;^c^ J dkWx^{LhdjWx^dy, by (I4.12p and (I4.22p . we obtain 

i,fc=l MiV 



N 



XI ^hcl J dkW^^{LhdjW^^)dy 



iV + 2 



j w^+^dyV^G{xo)cQ ■ Co + o{h^ 



(4.30) 
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For the integral / Ri{LhRi)dy, by K2T\) . Ri{ , ^ ) satisfies 



AR-R + pw^-'R + [V{xh) + A] ' 'm(a;;,)"f ^ dijV{xo)y^yjW 
-7:[Vixh) + \]* mixh)---'J2d,M^o)yryjwP = 0. (4.31) 



Hence 



R,{LhRi)dy = J R,iL,^R,)dy + Oih^) 
[l^(xft) + A] m{xh)~~ ^ 9ijl^(xo)9fc«V"(xo) / yiyjwLQ^{ykyiw)dy 

-4 iv ^ /" 

- [V'la;/,) + A] ^m(a;/,)"^"^ ^ 9ij\/(a;o)5fc/m(a;o) yiyjwLQ\ykyiwP)dy (4.32) 

1 _3 iv ^ /■ 

+ - [V^(x,i) + A] m{xh)~^~^ ^ dijm{xo)dkim{xo) yiyjwPLQ'^{ykyiwP)dy + 0{h^). 



i,j,k,l=l 

As in Section 3, we have used the following identities: 

AT 



^dijV{xo) j yiyjwljy=^ j |ypM;^^d?/Ay(xo) = 



^ dijV{xa) j y^yjW^^(l)2dy 

-[^(a;h) + A] ^m(x/,)"^ ^ 9ij1^(a;o)5fc/V'(xo) / yiyjwLQ^{ykyiw)dy 

1 _2 iv ^ /• 

- 2 t^*^^'"'* ^m(x/i)"^"^ ^ 9ijV"(xo)9fczm(xo) I yiyjwLQ^{ykyiw'P)dy , 



i,j,k,l=l 



f N 

E '9»ifc^(a;o)a:i,i / yjykwljy = ^ / |y|2w;2^ciyV(AV^)(a;o) ■ xi 

ij,fc=l IRAf JJJV 

TV 

E 9iife«'^^(a;o) / yiyjykyiwl^ = jv(j^+2) / li/l^'«^^;,(^i/^^'^(a^o) , 

i,j,k,l=l RJV Kiv 

TV 

E <9ijfcim(xo) / yiyjykyiwlt^ = n{n+2) I \y\^wPt^dyA^m{xo) 

i,j,k,l=l RiV RJV 
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RN KJV 

/ yl_iwLo\y%w)dy = ^ / r'^wL^^{r'^w)dy - ^^2(^+2) / r'^w^i{r)dy , 

R^ Rf^ R^ 

J yN^iyNwL^\yN^iyN'w)dy = jfj^:^ J r^w^i{r)dy , 

R!^ 

r / y%wL^\y%wndy = ^ / r^wL^^r^w^dy + J r'wMr)dy , 

RJV ]RiV ]Rjv 

/ yli^iwL^\yl,wP)dy = ^ / r^wL^\r^wP)dy - j^^j^^ J r^w^o{r)dy , 

M^f RJV KiV 

/ yN-iyNwLQ^{yN-iyNwP)dy = j^j^:^ J r^w^o{r)dy , 



where $0, '^'i satisfy 



"^0 + ^"^-o - "^-o + pwP-'% - ^$0 - r^wP = 0, r G (0, 00) 
<l>o(0) = $;,(0) = 0, 



and 



+ - $0 + - ^$1 -r'^w = 0,r e (0, cx)) , 

<l>i(0) = $UO) = 0, 

which can be proved as in Appendix B. 

Therefore, combining (giil, f09|) . fOOD and fO^jl . we obtain 

+ 0(1) = H,{xo) + Hsixo) + H,{xo) = H{xo) , (4.33) 

where 

g2(xo) = ^(^^^2) [n^o) + xr^ixo)-^ J \y\'w^d\y\)dy\\v'v{xo)\\l 

R'V 

^ [V{xo) + A] "'m(xo)^^-^ / \y\'w%{\y\)dyV^V{xo) ■ VM^o) 



N{N + 2) 
1 

4/V2 



+ 74^[^(^o) + A]^'m(xo)-^"^ / \y\'wPL^\\y\'wndy\Amixo)\' 



+ 2N{N + 2) [n^o) + A]"Mxo)^^-^ / |ypt.^<l>o(|2/|)dy||V^m(xo)|| 



2A^2(jv 



[l^(xo) + A] -'m{xo)-^-' J \y\'wP%{\y\)dy\Am{xo)\' , (4.34) 



28 



H3ixo)=^[Vixo) + \] 'm(xo)-^ j |?/| Wo??/V(Am)(xo) ■ xi 

J w^^'dyco-[V'G{xo)Y'co, (4.35) 



H^ixo) = g^^^ ^ [r(xo) + A] -^m(xo)-^ / WrfyAV(xo) 

^ [\/(xo) + A]~^m(xo)-^-W |2/|^tfP+MyA2m(xo). (4.36) 



8(iV + 2) 

Consequently, p{d") = 1 if H{xo) > 0, where H{xo) defined in (14.331) involves the i-th deriva- 
tives (for < z < 4) of \^ and m at xq. On the other hand, by Theorem 14.11 we have 
n{Lh) = n + 1. Thus we complete the proof of Theorem 11.41 by the orbital stability and 
instability criteria of [I7]-[18]. 

5 Appendix A 

In this section, we want to prove (I2.29P of Section 2, i.e. 

/i2 r 

{LhdjWccJ dkWx^dy = — — - / w^'j^^dydjkmixo) + o{h^) . (5.1) 

P + J- Jrn 

Proof. Note that by Lemma 12.31 and 12.41 we obtain 

LhdjW^^ = ^m{hy + Xh) - m{xh) pw^'^djiu^^ + m{hy + Xh)p{vl~^ - wl~^)djW^^ 

= — ^ diim{xo)yiyipwl~^djW^^ + h^m{xh)p{p - l)w^~V25jW^h + o(/i^) . 

i,l=l 

Hence we may write the integral Jj^jv {LhdjWx^) d^Wx^dy as follows: 

{LhdjW^J dkW^^ dy = Ii + l2 + o(/i^), (5.2) 



(5.3) 



where 

N 



h r 

h =ir dum{xo) / yiyipwH^djW^^dkW^Jy , 

^ ^,l=l -f^" 

h =h'^ / m{xh)p{p - l)wl~'^(j)2djW^^dkW^^dy . (5.4) 
Note that from (12. 3p . we have 

A - A + m{xh)pwll^ djkWx^ + m{xh)p{p - l)w^~^djW^^dkWr,^ = . (5.5) 
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Hence by f l2.14p . (15 ■4p and f l5.5p . we may use integration by parts to get 



h = -h' 



A — A + m{xh)pw^^ 



i-i 



2dy 



i,l=l 
'2 N 



^diim{xo) / ViViwP^djkW^^dy 



— y^diim{xo) I — ^dkW^^dy 



dum{xo) / yiyipwl^ ^djWxf^dkWxf^dy - h'^djkm{xo) / ykW^^dkHJ^f^dy 



i,l=l 
N 



- — <9i«m(xo) / yiyipwl^^djW^^dkW^^dy H —djkm{xQ) 



<:t'dy 



Combining (Q, (Q and (ES]), we obtain flSTT]) . 



6 Appendix B 

In this section, we prove f l3.27p . fl3.28p and fl3.29p of Section 3, i.e. 



j yWLo\y>ndy=^ J r'w^L-,\r'w^)dy 



2{N -1) 
iV2(iV + 2) 

yl^,w'''L-,\ylw^)dy = j r'w^L-,\r'w^)dy - ^ 



1 



r'^w^^o{r)dy , 



r'^w^^Q{r)dy , 



yN-WNW^LQ ^{yN~iyNwP)dy 



N{N + 2) 



r^w^^o{r)dy . 



where r := |?/| and $o satisfies 



$0 + ^^^"^0 - "^0 + pwP~^<^o - - r^njP = 0, r e (0, 



oo 



$o(0) = <f'o(0)=0. 
Proof. From (16. 4p . it is easy to check that 



ylfW^ , and Lq 



^ Vn-iVn 
— — 



yN-lVNW^ 



p 
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Then using the polar coordinate, we obtain 



yj^wPLf)\ylfwP)dy 



dy 



r cos 9n-iw^ 



^ . , r^cos^ 6'jv_i 1 ^ , X 1 r-i/ 2 «\ 
•^oW -2 -J^Mr) + j;^L,\r^wn 



dy 



J cos^ 6'Ar_i sin^ '^6j^_id6N-i 



/sin^-2^^_irf^^_i J 



TT 

/ COS^ 6'iv_i Sin^~^ 6'Ar„i(i6'Ar_i 

+ - — 

/sin^-2^^_^^^^_^ 



1 



r^w^^o{r)dy 



dy 



dy 



2.„P^-l..2.„.^.., , 2(iV-l) , , 



r w^^o{r)dy . 



This completes the proof of (16. ip . Similarly, one may obtain ( 16.2p and (]6.3p . respectively. 



7 Appendix C 

In this section, we prove f l4.12p of Section 4, i.e. 



J dkW^^Lh{djW^^ + htljj)dy 



N + 2 



wP+^dydjkG{xo)+o{h'). 
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Proof. Note that by (gS]), and ( KWf . we obtain 



and 



m{hy + Xh) - m{xh) pw^^ ^djW^^ 



+ m{hy + Xh)p{vl ^ - w^^ ^)djW^^ - V{hy + Xh) - V{xh) 



--h 
+ 



y ■ Vm{xh)pwl ^ + m{xh)p{p - - y ■ W{xh) 



Yl ^ii^^^hjViViPwl^ ^ + y ■ Vm(a;/,)p(p - l)wl^ Vi + rn{xh)p{p - V2 



TV 



1 1 V 

+ 7;^{xh)p{p - - 2X;;Vi - 9 XI 9ii^(^h)ym djw^^ + o(/i^) 



m{hy + Xh) - m(x/i) 



+ m(/i?/ + ^ - ty^^ - V{hy + x/,) - 



y ■ Vm{xh)pwl^ ^ + m{xh)p{p - ^4>i - y ■ "VVixh) 



+ h 



y ■ Vm{xh)pwl^ ^ + m{xh)p{p - - y ■ W (x/^) + 0{h^) 



Hence we may write the integral j dkWxf^LhidjWxf^ + hipj)dy as follows: 

J dkW^^Lh{djW^^ + hijj)dy = I0 + I1 + I2 + o(/i^), 



(7.2) 



where 



■ \/m{xh)pwl^ ^ + m{xh)p{p - l)w^^ - 1/ ■ VV^(x/i) ^jdkW^^dy , (7.3) 



„ 1 ^ 

Ji =/i^ y ^ [2 diim{xh)yiyipwll^ + ?/ ■ Vm{xh)p{p - l)wS;;Vi 



AT 



+ -m{xh)p{p -l){p- 2)wP^^(f)l - - X^i«^(3^ft)l/i?/« djW^^dkW^^dy , 



h=h' 



m{xh)p{p - l)tfS. ^hdjWx^dkWx^dy . 



(7.4) 
(7.5) 



Note that from ( 12. 3p . we have 

A - {V{xh) + a) + m{xh)pwl~^ djkWx,^ + m{xh)p{p - l)w^'^^djWx^dkWx^ = . (7.6) 
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Hence by f l4.6p . f l7.4p and (17. 5p . we may use integration by parts to get 



A - {V{xh) + A) + m{xh)pwl^ ^ djkW^Thdy 



A - {V{xh) + A) + m{xh)pwl^ ^ (f)2dy 



N 



y ■ W{xh)4>i - ^ X] 9iiV{xh)yiyiw^^ + y ■ Vm{xh)pw^^ Vi 



1,1=1 



N 



+ 9 ^ii^(^h)yiyiw'^^ + -m{xh)p{p - Vi djkW^^dy 



il=l 



N 



djV{xh)4>i + y ■ W{xh)dj(j)i + 7;^ diiV{xh)yiyidjW^^ + djkV{xh)ykWx^ 



1,1=1 



djm{xh)pwl^ Vi - y ■ Vm{xh)p{p - Vif^j^^x^ - 1/ ■ Vm{xh)pwl^ ^dj(t>i 



1 ^ 

- - ^ diim{xh)yiyipwll^djw^^ - djkm{xh)ykw 



p 



,1=1 



- -m{xh)p{p - - 2)wl^^(pldjW^^ - m{xh)p{p - Vi^i^i dkW^Jy 



h - 



y ■ Wm{xh)pwl^ + m{xh)p{p - <Pi-y- ^V{xh) djCpidkW^^dy 



djV{xh)(t)i + djkV{xh)ykWx^ - djm{xh)pwl Vi - djkm{xh)ykwl dkW^^dy 



(7.7) 



Note that from (14. 5p . we have 

A((9j0i) - [V'(xo) + A] d.j(t)i + m(xo)p<;;^9j>i + m{xo)p{p - l)wP~'^(j)idjW^g 
- y ■ W{xo)djW^o - djV{xo)w^^ + y ■ Vm{xQ)pwl~^djW^^ + djm{xo)wP^ = , 

and by direct computation, 



(7i 



LxoWxo = {p- ^)m{xo)wP^ , 
^xo(^w7xo + |y ■ Vtf^J = [V{xo) + A] Wx, 



xo 



(7.9) 



33 



Thus we may use f l7.3p - fl7.9p and integration by parts to get 
Io + h + h 



y ■ Vm{xh)pwl^ ^ + m{xh)p{p - l)w^^^(pi - y ■ W{xh) (i'j - dj(pi^dkW^^dy 
djV{xh)(pi + djkV{xh)ykWx^ - djm{xh)pwl~^(f)i - djkm{xh)ykW^^ dkW^^dy 



i)j - djcpi L,x^^pkdy + h / djV{xh) - djm{xh)pwl 



djkV{xh)ykWxh - djkm{xh)ykwl dkw^^^dy 



ndkW-x^dy 



djV{xo)wxg - djm{xo)wP^ ipkdy -h / djV{xh)w,j:,^ - djm{xh)w^^ dk(t)idy 



^djkV{xh)wl^ - --^djkm{xh)wP+^ dy + o{h^) 



djV{xo) {V{xo) + A) ^ (^yw^o + ■ Vw^J 



P 



^djkV{xQ)wl^ - ^^^9jfcm(xoX+^ dy + o{h?) 
djV{xQ){V{xQ) + A)~^(— + ]^y ■ Vw^^) 



P 



-mi^xo) ^djm{xo)wxo 9fcK(xo)tfxo - dkm{xo)w^g 



dy 



-djkV{xo)wl^ 



1 



djkm{xo)wP^^ dy + o{h 



p — 1 
1 1 



m(xo) '^djm{xo)dkV{xo) 



wl,dy 



-) [V{x^) + \]-^ d,V{xo)dkm{xo) 



p — 1 2p 

- ^^-^m(xo)^^9jm(xo)<9fcm(xo) / w^'^.^dy 



Recall that 



l:djkV{Xo)wl^ rl.,m('r-A^nP+'^ 

2 p 



-djkm{xo)w^'^^ dy + o(/i^) . 



(7.10) 



Wxo{y)= [V{xo) + A] ~ m{xoy^w{^/V{xo) + Xy) , 
m(xo) Vr(xo) = f [V{xo) + A] Vm(xo) , 

dijG{xo) = m{xo)~^~'^ m{xo)dijV{xo) + (1 - Y)9iV{xo)djm{xo) - f [V{xo) + A] dijm{xo) 
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and the integral identity 

Combining (Q and (ITrTOj) . we obtain ([71]). □ 
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